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ABSTRACT 

Repeated  surface  wave  group  velocity  measurements  over  the  same  paths 
were  made  by  using  large  aftershocks  of  several  great  (M  >  7)  earthquakes 
recently  occurred  in  China.  Multiple  filter  technique  was  applied  to  the 
properly  rotated  three-component  digital  data  from  Selsmological  Research 
Observatory  (SRO)  stations,  so  that  both  Rayleigh-  and  Love-wave  dispersion 
data  were  obtained  over  a  number  of  paths  crossing  various  tectonic  provinces 
of  China.  In  several  cases,  higher  mode  data  were  also  derived.  An  estimate 
of  uncertainty  of  these  dispersion  data  was  obtained  from  the  repeated  measure¬ 
ments  with  identical  source- receiver  geometry.  The  generalized  surface  wave 
Inversion  technique  was  applied  to  these  multi-mode  dispersion  data,  and 
crustal  and  upper  mantle  structures  were  derived  for  various  tectonic  provinces 
of  China.  The  results  clearly  demonstrate  that  the  Chinese  Mainland  is  far 
from  being  laterally  homogeneous,  the  lateral  heterogeneities  closely  reflect 

the  tectonic  developments  in  the  recent  past.  A  particularly  unusual  crystal 

c  r 

and  upper  mantle  structure  is  found  underlying  the  Tibet  plateau. 
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INTRODUCTION 

China  has  been  divided  tectonically  into  four  major  subplates  by  Sun 
and  Teng  (1977)  based  on  surface  geology.  These  subplates  include  (Figure  1): 

(1)  Chi nghai -Tibet  (southwestern  China), 

(2)  North  China, 

(3)  South  China,  and 

(4)  Northwestern  China  (including  the  Tien  Shan  fold  belt). 

These  subdivisions  are  similar  to  those  of  McElhinny  (1973),  based  on  paleo- 
magnetic  data;  Whittington  and  Hughes  (1972)  and  Jell  (1973),  based  on 
paleontologlc  data;  and  Dewey  and  Bird  (1970),  Hamilton  (1970)  and  Burrett 
(1974),  based  on  regional  geology.  Each  subplate  has  different  geomorphologic 
features,  gravity  anomalies  and  crustal  thickness  (Tung,  1974;  Sun  and  Teng, 
1977).  A  generalized  compilation  of  the  results  of  seismic  reflection  and 
refraction  as  well  as  gravity  surveys  has  resulted  in  the  inferred  crustal 
thickness  contours  shown  in  Figure  1  (Chinese  Academy  of  Sciences,  1974). 

Major  seismic  events  occur  alohg  the  subplate  boundaries.;  This  seismicity 
is  related  to  the  large-scale  tectonic  movements  discussed  by  Molnar  and 
Tapponler  (1975,  1977). 

The  division  of  China  into  four  major  subplates  is  helpful  from  the 
standpoint  of  surface  wave  dispersion  studies.  The  concept  of  regionalization, 
first  proposed  by  Tdksoz  and  Anderson  (1966),  consists  of  dividing  a  surface 
wave  path  Into  oceanic,  shield  and  mountain  tectonic  regions.  Each  region  has 
a  different  effect  upon  the  average  group  velocity.  Tbksoz  and  Anderson  (1966) 
noted  that  shield  areas  raise  the  average  group  velocity,  whereas  active 
tectonic  regions  have  the  effect  of  lowering  the  average  group  velocity.  As 
a  result,  surface  wave  paths  crossing  different  tectonic  regions  should  show 
different  trends  in  their  average  group  velocity  dispersion  curves.  The 
relationship  between  the  concept  of  regionalization  and  its  tectonic  implications 
has  been  studied  (Tung,  1974);  the  results  show  that  it  is  appropriate  to 


divide  China  tectonically  into  several  subplates,  so  that  relatively  pure 
path  surface  wave  dispersion  data  can  be  derived  for  crustal  and  upper 
mantle  structure  inversion.  This  paper  makes  use  of  the  digital  data  that 
has  become  available  from  numerous  large  aftershocks  due  to  the  recent 
occurrences  of  several  large  events  (M  >  7)  inside  China,  and  extracts  dis¬ 
persion  information  for  a  detailed  analysis  on  the  crustal  and  upper  mantle 
structures. 


Tectonic  Setting  of  China 

The  tectonic  setting  of  the  four  subplates  of  China  proposed  by  Sun 
and  Teng  (1977)  are: 

1)  The  Chinghai-Tibet  subplate  consists  of  the  area  north  of  the 
Himalayas  Including  the  Kun  Lun  Mountains.  Its  eastern  boundary 
coincides  with  the  western  border  of  Szechwan  and  Yunnan  provinces  at 

O 

approximately  longitude  105  E.  Its  western  boundary  is  less  well  de¬ 
fined,  extending  approximately  500  km  to  the  west  of .th^. western  border 
of  China.  The  borders  of  the  Chinghai-Tibet  subplate  are  coincident 
with  zopes  of  high  seismic  activity  and  fossil  subduction  (Molnar  and 
Tapponier,  1975,  1977). 

2)  The  North  China  subplate  is  bounded  to  the  north  by  the  Siberian 
platform  (Hamilton,  1970),  to  the  east  by  the  Pacific  Ocean,  to  the 
south  by  the  Tsinling  Shan,  and  to  the  west  by  the  Ordos  desert.  The 
Tslnling  Shan  consists  of  a  series  of  faults  striking  northwesterly  into 
the  Nanshan  fold  belt  and  Is  considered  as  both  a  distinct  paleontogical 
and  physical  boundary  between  the  two  subplates.  The  Shansi  graben  con¬ 
sists  of  a  series  of  northeast  striking  faults  which  extend  north  Into 
the  Ordos  platform.  These  features  are  evident  by  LANDSAT-1  imagery  as 
discussed  by  Cardwell  and  Isacks  (1976).  A  major  part  of  the  North  China 
subplate  is  the  North  China  plain,  a  sediment-filled  basin  with  an  average 


elevation  of  less  than  100  meters,  surrounded  by  uplifted  areas  of  more 
than  1  km  In  elevation.  Although  the  North  China  plain  is  considered  an 
"Intraplate"  area.  It  has  a  long  history  of  seismic  activity  as  Indicated 
by  the  Catalog  of  historical  earthquakes  In  China  (Lee,  et  al_-  1976). 

3)  The  South  China  subplate  is  south  of  the  Tsinling  Shan  and  east 
of  the  Tibetan  and  Tsaidam  plateaus.  It  consists  of  the  Yangtze  fault 
block  in  the  northwest  and  the  South  China  fold  block  in  the  southeast 
(Sun  and  Teng,  1977).  Its  eastern  boundary  runs  parallel  to  the  shelf 
margin  of  the  southeast  coast  of  China,  partially  coincident  with  the 
subduction  of  the  Pacific  plate  below  the  Eurasian  plate.  The  presence 
of  melange  in  the  east  Coast  Mountain  Range  of  Taiwan,  glaucophane  in 
the  Central  Mountain  Range  of  Taiwan  (Sun  and  Teng,  1977),  and  con¬ 
siderable  seismic  activity  in  northeastern  Taiwan  (Katsumata  and  Sykes, 
1969;  Tsai,  et  aT .  1978)  have  been  used  as  evidence  that  this  region  is 
an  active  zone  of  subduction. 

4)  The  Northwestern  China  subplate  consists  of  the-Tapim  and  Turfan 
basins,  the  Tien  Shan  fold  belt,  and  the  area  north  and  west  of  the  North 
China  subplate.  The  regional  tectonics  of  this  area  are  not  well  known. 

A  recent  interpreation  (Tapponier  and  Molnar,  1979)  based  on  Landsat 
imagery  and  seismicity  reveals  the  diversity  of  tectonics  in  the  western 
portion  of  this  subplate. 

Date 

The  earthquakes  used  in  the  present  study  (Table  1)  were  chosen  such  that 
a  major  portion  of  their  wave  paths  fall  within  the  four  major  subplates  of 
China  (Fig.  2).  A  number  of  earthquakes  are  used  for  each  path  in  order  to 
obtain  an  estimate  on  the  data  repeatability  and  observational  uncertainties. 
The  path  lengths  range  from  approximately  1700  km  to  over  5000  km.  For  the 
path  Tangshan  to  Taipei,  the  accuracy  of  surface  wave  data  for  periods  greater 
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than  64  seconds  was  limited.  All  events  in  this  study  were  recorded  at 
two  Seismological  Research  Observatory  (SRO)  stations  located  in  Mashad, 

Iran  (MAIO)  and  Taipei,  Taiwan  (TATO). 

SRO  data  in  digital  form  on  tape  is  convenient  for  computer  processing. 

The  SRO  stations  operate  at  high  sensitivity  levels,  thus  expanding  the 
detection  threshold  for  seismic  events  and  making  small -magnitude  events 
(M  -v  5)  useful  for  long-  and  short-  period  surface  wave  dispersion  analysis. 
Each  seismic  event  must  be  decoded  and  plotted  from  the  data  tapes.  Pure 
Love  and  Rayleigh  waves  are  obtained  by  rotating  the  data  coordinates.  An 
example  of  a  plotted  seismic  event,  unrotated  and  rotated,  is  given  in  Figs. 

3  and  4.  A  compllete  discussion  of  handling  SRO  data  can  be  found  in  Peterson, 
et  al_.  (1976). 

Method 

A  multiple  filtering  analysis  (Dziewonski  and  Hales,  1970;  Herrmann,  1973; 
Tung,  1974;  Seekins  and  Teng,  1976)  has  been  applied  to  the  surface-wave  train 
to  obtain  group  velocity  dispersion  curves.  A  typical  contoilred  result  of  the 
relative  amplitude  of  wave  energy  arrivals  over  the  path  Tangshan  to  Taipei 
is  shown  for  both  the  Rayleigh  and  the  Love  components  in  Figs.  5,  6  and  7. 

In  these  figures  the  dots  represent  the  average  group  velocities  obtained 
from  the  multiple  filtering  analysis.  At  the  low-period  end  (below  20  seconds) 
ambiguity  often  arises  in  drawing  contours  due  to  interference  by  S  waves. 

The  contour  lines  will  give  several  local  maxima  and  will  not  allow  a  clearcut 
determination  of  group  velocity  values.  The  average  dispersion  curves  for 
Rayleigh  and  Love  components  are  plotted  for  all  events  of  each  path  in 
Fig.  8-15.  Higher-mode  data  are  included  when  they  are  observable  from  the 
multifilter  analysis.  The  observed  data  over  each  path  are  listed  in  Tables 
2,  3,  4  and  5. 
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The  Inversion  of  the  surface  wave  dispersion  data  for  determining 
earth  structure  is  based  on  the  Backus-Gilbert  inversion  method  (1967, 

1968,  1970).  Oetailed  descriptions  of  this  method  have  been  presented 
(Jackson,  1972;  Wiggins,  1972;  and  Crosson,  1976).  The  basic  logic  of 
the  inversion  theory  and  some  experience  in  its  numerical  applications 
will  be  briefly  discussed  as  follows. 

The  observations,  ,  are  related  to  certain  model  paramters,  X^, 
in  some  known  way,  Y^  *  A(Xi,  X21  .  .  .  X  ).  A  quasi-linear  relationship 
is  established  by  taking  a  Taylor  series  expansion  about  some  initial  value 
X'  . 


Y*  a  (X1)  +  i  AXi  +  higher  order  terms 

1  1  0  ”3X7  >X'  J 

J  • 

* 

where  Y.  *  A-  (X*)  +  AY. 

lip  I 

Ignoring  the  higher  order  terms,  we  obtain  the  function  relationship  in 
matrix  form: 

AAX  «  AY 

where  AY^  »  observed  data  minus  calculated  values,  and  A  is  a  matrix  whose 

elements  are  the  partial  derivatives  in  Equation  1.  AX  is  a  p  x  1  vector 

whose  elements  are  corrections  to  the  initial  model  X'.  AY  is  the  difference 

o 

between  the  observed  and  theoretical  values. 
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In  the  present  case  the  observations  are  group  velocities  of  Love  and  Rayleigh 
waves,  and  the  model  parameters  are  shear  velocities  and  densities  for  a  horizon¬ 
tally  layered  earth.  The  partial  derivatives  in  equation  1  are  generated  numeri¬ 
cally  by  the  method  described  by  Rodi  et  al_.  (1975). 

The  classic  least  squares  solution  to  equation  2  is  by  minimizing  the 
Euclidean  length  of  A AX -AY,  or 

1|  AAX-AY  ||  =  (AAX-AY)T(AAX-AY)=  e 

Given  equation  2,  the  solution  is  (Hanson  and  '  ;wson,  1975): 

AAX  =  AY 
ATAAX  =  ATAY 
AA  =  (AVlATAY 

where  AX  is  the  estimated  correction  vector. 

A  suitable  inverse  to  this  problem  can  be  obtained  by  using  ?the  singular 

f 

value  decomposition  (SVD)  of  the  matrix  A. 

A  =  UAVT 

where  the  columns  of  U  are  the  eigenvectors  associated  with  the  columns  of  A, 
the  rows  of  V  are  the  eigenvectors  associated  with  the  rows  of  A,  and  A  is  a 
diagonal  matrix  of  non-zero  eigenvalues  of  the  matrix  A. 
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If  we  let 

H  *  <aW 


(3a) 


then 

HAM  *  HAY  (3b) 

and 

AX  =  HA  X  (3c) 

and 

AX  =  HAY 

Using  the  SVD  we  obtain  (Jackson,  1972): 

AX  =  (AV'V)AY  (3d) 

H  *  AV'V 

*  f  r 

Here,  H  is  known  as  the  pseudo-inverse  of  equation  2. 

The  matrix  product  HA=R  is  known  as  the  resolution  matrix,  and  is  a  measure 
of  the  uniqueness  of  the  solution  (Jackson,  1972).  The  rows  of  R  are  called  the 

A  A 

resolving  kernels.  From  equation  3c  it  can  be  seen  that  an  element  of  AX,  AX^, 
may  be  interpreted  as  convolving  the  i^  row  of  R  with  the  vector  AX.  Therefore, 
M.  can  be  considered  the  weighted  sum  of  nearby  values  (Jackson,  1972). 

Because  ATA  is  nearly  singular,  a  number  of  problems  arise:  (1)  The  solution 
vector  becomes  large,  which  can  cause  the  problem  to  leave  the  region  of  linearity; 
(2)  the  solution  oscillates  with  each  iteration.  By  looking  at  the  pseudo-inverse 
H  we  can  see  the  cause  of  this  instability. 
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H  -AvV 

_  1 

where  the  matrix  A  equals 


_  i 

A  * 


and  the  A.'s  are  the  eigenvalues  of  the  matrix  ATA.  If  ATa  is  near  singular, 
one  or  more  of  the  eigenvalues  will  be  approaching  zero.  From  equation  3  it  can 
be  seen  that  a  small  eigenvalue  will  cause  a  large  change  in  one  or  more  values 
of  the  correction  vector  AX.  Similarly,  it  has  been  shown  that  the  variance  of 

*•  •  f  • 

the  model  parameters  are  inversely  proportional  to  A  (Jackson,  1972).  Thus,  small 
A's  produce  large  standard  deviations. 

There  are  two  main  approaches  to  stabilizing  the  inversion  process.  One 
method  is  to  examine  the  eigenvalue  spectrum  of  the  matrix  A.  Small  eigenvalues 
are  removed  when  the  variance  becomes  too  large,  but  removing  the  eigenvalues 
degrades  the  resolution.  This  trade-off  between  resolution  and  variance  has  been 
discussed  by  a  number  of  authors  (Jackson,  1972;  Wiggins,  1972;  and  Crosson,  1976). 
Braille  and  Keller  (1975)  used  this  method  for  the  inversion  of  group  velocity 
data. 

Inversion  stability  can  also  be  achieved  by  using  Marquardt's  method 
(Harquardt,  1963)  or  the  stochastic  inverse  (Franklin,  1970).  These  procedures 
lead  to  the  suppression  of  small  eigenvalues.  Franklin's  stochastic  inverse 
estimates  the  parameter  in  the  presence  of  noise. 


AAX  +  n  =  AY  (4) 

Where  n  is  a  vector  of  observational  noise,  Franklin  (1970)  shows  that  the 
solution  to  equation  4  is  given  by: 

AX  *  WAT(AWAT  +  EOJ^AY  (5) 

where  W  is  the  covariance  matrix  of  the  parameter  and  ED  is  the  covariance  matrix 
of  the  observations. 

Marquardt's  method  considers  a  minimization  of  the  functional 
(AAX  -  AY)T(EDT£D)  (AAX  -  AY)  -  AxVw)  AX  (6) 

with  the  solution  given  by: 

ax  =  (ated_1  a  +  w  ‘Wed"1  y  (7) 

T  _  i  T 

Letting  A  ED  =  A  ,  equation  7  becomes: 

AX  =  (ATA  +  o2I)  ATAY  ‘  {  f  (8) 

2 

where  a  is  the  variance  of  the  parameter. 

By  looking  at  the  $VD  of  equation  8  we  can  see  how  Marquardt’s  method  tapers 
the  eigenvalue  spectrum: 

AX  *  V|  (A*  +  o2I)"1A|UTAY  (9) 

An  element  in  brackets  from  equation  9  equals:  (Crosson,  1976). 

|(A2  +  02I)"lA|.  =  A./(A2.  +  a) 

* 

So  as  X j  approaches  zero,  the  parameter  AX^  also  approaches  zero.  A  similar 
analysis  for  the  stochastic  inverse  shows  that  it  too  suppresses  small  eigenvalues. 


The  quantity  a2  can  also  be  thought  of  as  a  tradeoff  parameter  between 
resolution  and  variance  (Der  et  aK .  1970;  Wiggins,  1972;  Crosson,  1976).  In 
Marquardt's  method,  the  resolution  equals  (Crosson,  1976): 

R  «  HA  *  (ATA  +  a2I)ATA  -  V|  (A  +  x^PA^V1  (10) 

For  the  stochastic  inverse 

R  -  |  (A2  +  o2I)"*|A2WT  (ID 

The  quantity  in  brackets  of  equations  10  and  11  is  a  diagonal  matrix  whose 
elements  have  the  form: 

Xi2/(X2  +  o2) 

If  a2  *  0,  then  R  =  I,  but  as  o2  approaches  zero,  then  the  variance  of  the 
parameters  becomes  larger.  We  then  adjust  the  size  of  o2  until  there  is  an 
acceptable  trade-off  between  standard  deviation  and  resolution. 

Two  other  elements  that  must  be  analyzed  in  the  i version  process  are  the 
apriori  (SDX0)  and  a  posteriori  (SDX)  standard  deviations  of  the  inversion 
parameters: 

SDX0  = 

SDX  = 

where  D  =  a  row  vector  of  a  delta  matrix  a  (0,0,. . . ,1 ,0,.. . ,0).  This  analysis 
will  give  a  measure  of  how  much  new  information  is  coming  from  the  data.  SDX0 
is  the  standard  deviation  of  the  initial  model  paramters  before  inversion, 
whereas  SDX  is  the  standard  deviation  calculated  from  the  data  after  inversion. 
If  SDX0  and  SDX  are  approximately  the  same,  we  know  that  very  little  additional 
Information  is  coming  from  the  data. 

The  advantages  of  Marquardt's  method  and  the  stochastic  inverse  are  that 
(1)  the  SVD  of  A  is  not  explicitly  determined,  and  (2)  there  is  no  decision  to 
be  made  on  the  range  of  the  matrix  A^A. 
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A  combination  of  these  two  methods  is  used  in  our  inversion  process.  When 
the  number  of  parameters  Is  greater  than  the  number  of  data,  the  stochastic 
Inverse  Is  used.  If  the  number  of  observations  is  greater  than  the  number  of 
parameters,  Marquardt's  method  Is  used.  This  procedure  leads  to  the  minimum 
number  of  computations. 

Results  and  Discussion 

Initial  models  for  the  Inversion  were  based  on  the  results  of  a  previous 
surface  wave  study  done  by  Tung  (1974).  Starting  models  for  all  paths  were 
initially  inverted  for  both  shear  velocity  and  density.  During  the  inversion, 
large  instabilities  of  the  density  parameter  were  found  in  the  crust.  This  was 
corrected  by  changing  the  starting  model  and  inverting  for  the  compressional 
wave  velocity  with  the  shear  velocity  and  density.  As  discussed  by  Wiggens 
(1972),  Oer,  et  aK  (1978),  and  Jackson  (1976),  the  partial  derivatives  for 
shear  velocity  are  larger  than  those  for  compressional  velocity  or  density. 
Therefore  the  resolution  for  shear  wave  velocity  is  considerably  better  than 

the  other  two  parameters.  Only  the  resolution  matrix  for  the  partial  de- 

;  f 

rivatives  with  respect  to  shear  velocity  shall  be  presented. 

1 .  Chi nghai -Tibetan  Plateau  and  Hindukush  Region 

Two  separate  paths  across  the  Tibetan  plateau  were  used  (Fig.  2).  Both 
paths  originate  from  earthquakes  located  on  the  eastern  boundary  of  the  Tibetan 
plateau.  Approximately  3/4  of  the  paths  encompass  the  Chi nghai -Tibetan  subplate. 
The  first  group  occurred  in  Szechwan  province  and  the  second  group  occurred 
to  the  south  in  Yunnan  Province.  A  list  of  the  surface  waves  is  presented 
in  Table  1.  For  each  path  across  the  Tibetan  plateau  inversions  were  performed 
using  the  fundamental  and  first  higher  mode  Rayleigh  waves.  A  fit  of  the 
theoretical  model  to  the  observed  group  velocity  is  presented  in  Figs.  16  and  17. 

A  comparison  between  the  average  continental  dispersion  for  Rayleigh  waves  (Ewing, 
et  al_.  1957)  and  Rayleigh  waves  for  the  four  paths  across  China  are  presented  in 
Fig.  18.  The  group  velocities  across  the  Tibetan  plateau  are  substantially  lower 
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than  the  group  velocities  for  the  average  continental  structure,  with  differences 
extending  to  periods  of  120  seconds,  indicating  abnormal  upper  mantle  structure. 

Group  velocities  across  the  Tibetan  plateau  are  the  lowest  of  all  the 
regions  of  China.  Rayleigh  wave  group  velocities  ranged  between  2.60  to  3.65  km/sec 
and  steep  slopes  occurred  in  the  dispersion  curves  between  40  to  60  seconds,  where 
the  change  in  group  velocity  was  as  high  as  0.6  km/sec.  Love  waves  were  not 
used  in  the  inversions  due  to  the  absence  of  amplitude  maxima  in  the  dispersion 
curves  between  30  to  60  seconds.  Earlier  studies  (Thatcher  and  Brune,  1969;  James, 

1971)  have  discussed  the  problem  of  interference  of  Love  waves  crossing  inhomogeneous 
media.  Tung  (1974)  attributed  this  problem  on  the  Tibetan  plateau  to  Love  waves 
following  a  non-least  time  path.  A  possible  structural  implication  is  that  there 
may  be  structures  along  the  path  corresponding  to  wave  lengths  of  the  Love  waves, 
causing  the  wave  energy  to  travel  along  a  non-least  time  path  with  slower  velocities. 

Starting  and  final  shear  wave  velocity  models  with  standard  deviation  bars 
for  the  path  Szechwan  to  Mashad  are  presented  in  Figure  19.  Table  6  presents 
the  initial  and  final  model  parameters  and  the  resolution  matrix^for  each  model. 

The  resolution  matrix  for  the  eight-layer  model  is  presented  in  Table  7.  From 
the  width  of  the  resolving  kernels  it  is  possible  to  resolve  layers  of  15  km  in 
the  upper  crust,  15-20  km  in  the  middle  crust,  and  25-30  km  in  the  lower  crust. 

Starting  and  final  shear  wave  velocity  models  with  standard  deviation  bars 
for  the  path  Yunnan  to  Mashad  are  presented  in  Figure  20.  Tables  8-9  present 
the  initial  and  final  model  parameters  and  the  resolution  matrix.  From  the  width 
of  the  resolving  kernels  it  is  possible  to  resolve  a  10-15  km  layer  at  the  top 
of  the  crust,  15-20  km  layer  in  the  middle  of  the  crust,  and  a  20-25  km  layer 
at  the  base  of  the  crust. 

Resolution  in  the  mantle  is  poor  for  both  paths.  The  standard  deviations 
were  found  to  be  less  than  0.15  km/sec  for  all  the  models  presented.  Increasing  ! 

the  resolution  in  the  crust  and  mantle  leads  to  both  large  standard  deviations 
and  Instabilities  in  the  shear  wave  velocity  model  (Rosenthal,  1977).  A  measure 

S 
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of  the  amount  of  new  information  obtained  from  the  data  is  given  by  comparing 
the  magnitudes  of  the  apriori  standard  deviation  (SDX0)  to  the  a  posteriori 
standard  deviations  (SOX)  given  in  Tables  7  and  9.  In  the  upper  mantle  there 
is  only  a  small  amount  of  new  Information  being  obtained  from  the  data,  and 
the  final  model,  including  the  location  of  the  Moho,  is  most  affected  by  the 
starting  model.  As  a  result,  a  low  velocity  layer  in  the  upper  mantle  can 
only  be  inferred  from  the  small  amount  of  new  information  being  obtained. 

Additional  data  for  periods  longer  than  120  seconds  and  good  Love  wave  data 
would  significantly  increase  the  amount  of  new  information. 

The  shear  wave  velocity  structure  for  both  paths  crossing  the  Tibetan 
plateau  is  similar.  The  crust  is  approximately  70  km  thick  and  may  be  broken 
into  four  layers.  The  first  layer  has  shear  velocities  that  range  between 
2.8 7-3.06  km/sec,  suggesting  a  layer  approximately  15  km  thick  of  partly  sedi¬ 
mentary  origin.  Tung  (1974)  and  Mu,  et  al^.  (1974)  discussed  the  existence  of 
thick  sedimentary  layers  in  the  Mount  Jolmo  Lungma  (Mount  Everest)  region  of 
the  Himalayas.  Both  Birch  (19§3)  and  Simmons  (1964)  reported  shear  velocities 

i  f  ' 

for  sedimentary  rocks  at  10  to  15  km  depth  of  approximately  3.00  to  3.08  km/sec. 

A  second  layer  has  shear  wave  velocities  ranging  from  3.20  to  3.40  km/sec  and 
is  approximately  15  km  thick,  suggesting  a  layer  of  granitic  or  equivalent  com¬ 
position.  A  third  layer  has  shear  velocities  between  3.40  to  3.70  km/sec.  Is 
approximately  15  to  20  km  thick,  and  also  seems  to  consist  of  granite.  For 
comparison,  Birch  and  Bancroft  (1958),  Birch  (1963)  and  Simmons  (1964)  reported 
shear  velocities  for  granite  under  3  to  5  kilobars  of  pressure  of  3.45-3.70  km/sec. 
Simmons  (1964)  reported  that  granite  with  a  shear  wave  velocity  of  3.56  km/sec 
at  4  kilobars  pressure  has  a  shear  wave  velocity  of  3.79  km/sec  at  pressures  of 
10  kilobars.  The  fourth  layer  with  shear  wave  velocities  between  3.85  to  3.95 
km/sec  corresponds  to  granite  at  pressures  of  approximately  20  kilobars,  which 
would  correspond  to  depths  of  60  to  70  km. 

The  dispersion  curves  and  shear  wave  velocity  profile  for  the  Tibetan 
plateau  is  in  agreement  with  those  values  presented  by  Tung  (1974)  and  Bird 
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of  the  amount  of  new  information  obtained  from  the  data  is  given  by  comparing 
the  magnitudes  of  the  apriori  standard  deviation  (SDX0)  to  the  a  posteriori 
standard  deviations  (SDX)  given  in  Tables  7  and  9.  In  the  upper  mantle  there 
is  only  a  small  amount  of  new  Information  being  obtained  from  the  data,  and 
the  final  model,  including  the  location  of  the  Moho,  Is  most  affected  by  the 
starting  model.  As  a  result,  a  low  velocity  layer  in  the  upper  mantle  can 
only  be  inferred  from  the  small  amount  of  new  information  being  obtained. 

Additional  data  for  periods  longer  than  120  seconds  and  good  Love  wave  data 
would  significantly  increase  the  amount  of  new  information. 

The  shear  wave  velocity  structure  for  both  paths  crossing  the  Tibetan 
plateau  Is  similar.  The  crust  is  approximately  70  km  thick  and  may  be  broken 
into  four  layers.  The  first  layer  has  shear  velocities  that  range  between 
2.87-3.06  km/sec,  suggesting  a  layer  approximately  15  km  thick  of  partly  sedi¬ 
mentary  origin.  Tung  (1974)  and  Mu,  et  al_.  (1974)  discussed  the  existence  of 
thick  sedimentary  layers  in  the  Mount  Jolmo  Lungma  (Mount  Everest)  region  of 

the  Himalayas.  Both  Birch  (19$3)  and  Simmons  (1964)  reported  shear  velocities 

if 

for  sedimentary  rocks  at  10  to  15  km  depth  of  approximately  3.00  to  3.08  km/sec. 

A  second  layer  has  shear  wave  velocities  ranging  from  3.20  to  3.40  km/sec  and 
is  approximately  15  km  thick,  suggesting  a  layer  of  granitic  or  equivalent  com¬ 
position.  A  third  layer  has  shear  velocities  between  3.40  to  3.70  km/sec,  is 
approximately  15  to  20  km  thick,  and  also  seems  to  consist  of  granite.  For 
comparison,  Birch  and  Bancroft  (1958),  Birch  (1963)  and  Simmons  (1964)  reported 
shear  velocities  for  granite  under  3  to  5  kilobars  of  pressure  of  3.45-3.70  km/sec. 
Simmons  (1964)  reported  that  granite  with  a  shear  wave  velocity  of  3.56  km/sec 
at  4  kilobars  pressure  has  a  shear  wave  velocity  of  3.79  km/sec  at  pressures  of 
10  kilobars.  The  fourth  layer  with  shear  wave  velocities  between  3.85  to  3.95 
km/sec  corresponds  to  granite  at  pressures  of  approximately  20  kilobars,  which 
would  correspond  to  depths  of  60  to  70  km. 

The  dispersion  curves  and  shear  wave  velocity  profile  for  the  Tibetan 
plateau  is  In  agreement  with  those  values  presented  by  Tung  (1974)  and  Bird 
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and  Tbksoz  (1977).  Chun  and  Yoshl  (1977),  using  Rayleigh  and  Love  waves 
for  many  different  paths  across  the  Tibetan  plateau,  reported  a  crustal  thick¬ 
ness  of  approximately  70  km  and  observed  a  velocity  reversal  in  the  middle 
crust.  Due  to  a  15  to  20  km  resolution  in  the  middle  crust,  it  is  doubtful 
a  velocity  reversal  could  be  observed  for  data  over  the  two  paths  presented. 

Due  to  the  small  ampunt  of  information  (apriori  vs.  a  posteriori  standard 
deviations)  derived  from  inverting  Rayleigh  waves  in  the  upper  mantle,  the 
above  structural  conclusions  are  only  tentative.  The  upper  mantle  shear 
velocities  are  lower  than  either  the  Gutenburg  earth  model  (Tacheuchi,  et  al. 

1964)  or  the  Brune  and  Dorman  (1963)  model  for  the  Canadian  Shield  (Fig.  21). 

The  velocity  at  the  top  of  the  mantle  is  approximately  4.4  km/sec,  with  in¬ 
dications  of  a  low  velocity  layer  at  about  90  km  depth.  Since  the  low  velocity 
layer  can  only  be  inferred,  the  existence  of  extremely  low  velocities  of  4.1 
km/sec  is  tenative.  Low  shear  velocities  in  the  upper  mantle  of  Tibet  have 
been  reported  by  Tung  (1974),  which  suggests  a  high  degree  of  partial  melting 

and  a  low  resistance  to  deformation  (Anderson,  et  al_.  1972). 

»  f 

2.  North  and  Northwestern  China 

The  path  between  Tangshan  and  Mashad,  Iran  (Fig.  2)  crosses  the  North 
China  subplate  and  the  Northwestern  China  subplate  including  the  Tien  Shan 
fold  belt.  Tung  (1974)  found  crustal  thicknesses  of  30,  40  and  45  km  for 
Northeast,  Central  and  Northwest  China  respectively.  Crustal  thicknesses  in 
the  Tien  Shan  and  Pamir  Mountains  have  been  reported  as  high  as  65  km  (Arkhangel’ 
skaya,  1964). 

A  list  of  the  earthquakes  and  surface  wave  trains  used  over  the  path  Tangshan 
to  Mashad,  Iran  are  listed  in  Table  1.  The  fit  of  the  theoretical  to  observed 
group  velocities  for  the  Rayleigh  and  Love  waves  is  given  in  Fig.  22.  The 
theoretical  Rayleigh  wave  group  velocity  is  lower  than  the  average  continental 
Rayleigh  wave  dispersion  (Ewing,  e_t  al_.  1957).  Starting  and  final  shear  wave 
velocity  models  with  standard  deviation  bars  are  presented  in  Fig.  23.  Tables 
10  and  11  present  the  starting  and  final  model  parameters  and  the  resolution 


16 


matrix.  From  the  width  of  the  resolving  kernels  presented  in  Table  11,  it  is 
possible  to  resolve  a  10  km  layer  at  the  top  of  the  crust,  10  to  15  km  layer 
in  the  middle  of  the  crust,  and  a  20  km  layer  at  the  base  of  the  crust. 

When  Rayleigh  and  Love  waves  are  inverted  simultaneously,  a  considerable 
amount  of  new  Information  may  be  derived  from  the  data.  Observing  the  magni¬ 
tude  of  the  apriori  (SDX0)  vs.  the  a  posteriori  (SDX)  standard  deviations,  a  con¬ 
siderable  amount  of  new  information  Is  obtained  at  depths  of  100  km  (seventh 
layer).  Below  this  depth,  less  new  information  is  obtained  from  the  data,  and 
the  final  model  is  most  affected  by  the  starting  model. 

The  average  crustal  thickness  over  the  path  from  Tangshan  to  Mashad  is 
approximately  45  km.  Since  this  path  is  not  considered  a  pure  path,  only  average 
shear  velocities  are  presented  (Fig.  2).  A  three-layer  crust  seems  to  be  an 
adequate  model  for  this  path.  The  first  layer  is  approximately  10  km  thick  and 
has  an  average  shear  velocity  of  approximately  3.0  km/sec,  suggesting  a  layer 
at  least  partly  sedimentary  (Birch,  1963;  and  Simmons,  1964).  The  second  layer 
is  15  to  20  km  thick  and  has  arv  average  shear  wave  velocity  <jf  3.55  km/se"c.  The 
layer  at  the  base  of  the  crust,  20  to  25  km  thick,  has  an  average  shear  velocity 
of  approximately  3.80  km/sec,  suggesting  a  layer  of  granitic  origin  (Birch,  1963; 
Simmons,  1964).  The  average  shear  wave  velocity  at  the  top  of  the  mantle  is  4.40 
to  4.45  km/sec  and  the  low  velocity  layer  begins  at  approximately  90  km.  Using 
a  three-layer  earth  model  proposed  by  Dorman  (1959),  Schechkov  (1961)  was  able  to 
fit  mixed  path  dispersion  data  for  Northwest  China  to  a  theoretical  three-layer 
model  of  average  crustal  thickness  of  50  km;  the  first  layer  has  a  shear  velocity 
of  3.40  km/sec  and  is  20  km  thick,  the  second  layer  has  a  3.53  km/sec  shear  velocity 
and  is  30  km  thick,  and  a  semi -infinite  third  layer  has  a  4.50  km/sec  shear  velocity 
for  the  mantle.  Applying  the  concept  of  regionalization  to  paths  crossing  North¬ 
west  China,  Tung  (1974),  presented  a  three-layer  crust  of  45  km  thickness.  Shear  I 

i 

velocities  calculated  by  Tung  (1974)  for  the  crust  and  upper  mantle  are  lower  than  jj 
the  mixed  path  average  shear  velocities  presented.  The  low  velocity  zone  presented  f 
is  not  as  prominent  as  Tung  (1974)  reported,  and  can  only  be  inferred  from  the 
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information  obtained  from  the  data. 

A  comparison  of  the  Canadian  Shield  model  (Brune  and  Dorman,  1963)  with 
the  average  shear  velocity  model  for  Tangshan  to  Mashad  is  presented  in  Fig. 

24.  Both  models  have  higher  shear  velocities  in  the  crust  and  at  the  top  of 
the  mantle  than  the  theoretical  model  presented  for  the  path  Tangshan-Mashad. 

The  tectonics  over  the  path  between  Tangshan-Mashad  may  account  for  the  lower 
velocity  at  the  top  of  the  mantle.  The  catalog  of  historical  earthquakes  pre¬ 
sented  by  the  Chinese  Institute  of  Geophysics  Academia  Sinica  (Lee,  et  al^.  1976) 
shows  considerable  seismic  activity  throughout  the  regions  traversed  by  the 
surface  waves.  Successive  stages  of  subduction  of  the  Indian  plate  below  the 
Eurasian  plate  are  believed  to  have  taken  place  in  areas  as  far  north  as  the 
Tarim  basin  and  Mongolia  (Chang  and  Zeng,  1973;  Zonenshain,  1973).  These  are 
consistent  with  a  thicker  crust,  lower  group  velocities  and  lower  shear  velocities 
in  the  upper  mantle.  These  results  suggest  that  the  origin  of  the  Tien  Shan  fold 
belt  may  also  be  a  result  of  the  collision  of  the  Eurasion  and  Indian  plates. 

3.  Coastal  China  %  *  f 

For  most  of  the  path  between  Tangshan  to  Taipei  the  surface  waves  are 
parallel  to  the  east  coast  of  China.  According  to  the  tectonic  map  of  the 
Institute  of  Geology  of  The  Chinese  Academy  of  Sciences  from  which  Fig.  1  is 
derived,  the  crustal  thickness  along  most  of  the  path  is  about  30  km.  A  starting 
model  for  inversion  was  chosen  to  reflect  the  30  km  crustal  thickness.  A  list 
of  the  earthquakes  and  surface  wave  trains  used  in  the  analysis  is  given  in  Table 
1.  The  fit  of  the  theoretical  to  the  observed  Rayleigh  and  Love  wave  group 
velocities  is  presented  in  Fig.  25.  A  comparison  between  the  theoretical  Rayleigh 
and  the  average  continental  Rayleigh  wave  group  velocities  (Ewing,  et^al_.  1957) 
shows  lower  group  velocities  of  both  short  and  intermediate  periods. 

The  starting  and  final  shear  wave  velocity  models  with  standard  deviation 
bars  are  presented  in  Fig. 26.  The  shear  velocity  parameters  and  resolution  matrix 


are  presented  in  Tables  12  and  13.  From  the  width  of  the  kernels  of  the 
resolution  matrix  (Table  13),  it  is  possible  to  resolve  a  5  to  7  km  layer  at 
the  top  of  the  crust,  a  10  km  layer  in  the  middle  of  the  crust,  and  a  15  to 
20  km  layer  at  the  base  of  the  crust.  An  examination  of  the  apriori  (SDX0) 
and  a  posteriori  (SOX)  standard  deviations  indicates  that  a  small  amount  of 
information  exists  in  the  data  below  the  fifth  layer.  Thus  longer  period  data 
is  needed  for  increased  information  in  the  mantle.  Since  the  path  length 
from  Tangshan  to  Taipei  is  approximately  1700  km,  with  the  present  analysis 
periods  greater  than  64  seconds  will  not  have  proper  mode  separation.  Con¬ 
sequently,  it  is  unlikely  that  for  this  path  new  information  below  the  top  of 
the  mantle  can  be  obtained  and  shear  velocity  in  the  mantle  is  most  affected  by 
the  starting  model. 

A  crust  with  a  thickness  of  30  km  and  composed  of  three  layers  ia  an  adequate 

model  for  the  structure  over  the  path.  The  first  layer,  approximately  5  to  7  km 

thick,  seems  to  be  a  sedimentary  layer  with  shear  velocity  of  approximately  2.95 

km/sec.  The  second  layer  is  'v  JO  km  thick.  The  third  layer,  approximately  15 
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to  20  km  thick,  corresponds  to  s  shear  velocity  of  approximately  3.85  km/sec.  in 
the  mantle,  a  low  velocity  layer  can  only  be  inferred  due  to  the  lack  of 
sufficient  information  from  the  data.  A  comparison  of  the  Canadian  Shield  model 
(Brune  and  Dorman,  1963)  and  the  Gutenberg  earth-model  (Tacheuchi ,  et  al_.  1964) 
with  the  shear  velocity  profile  shows  higher  velocities  in  the  mantle  and  the 
base  of  the  crust  (Fig.  27). 

The  shear  velocity  of  4.46  km/sec.  at  the  top  of  the  mantle  is  lower  than 
the  shear  velocity  of  the  Canadian  Shield  (4.65  km/sec.).  Tung  (1974)  also 
calculated  a  mantle  shear  velocity  of  4.42  km/sec.  at  the  top  of  the  mantle  for 
a  path  crossing  Northeast  China  recorded  in  Hong  Kong.  The  tectonics  over  the 
path  from  Tangshan  to  Taipei  explains  the  lower  velocities  in  the  mantle. 

The  region  traversed  is  considered  an  "intraplate"  region,  although  a  con¬ 
siderable  amount  of  seismic  history  is  evident  from  the  catalog  of  historical  earth¬ 
quakes  compiled  by  the  Chinese  Institute  of  Geology  Academia  Sinica  (Lee,  et  al.  197 
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For  example,  in  the  northeastern  portion  of  the  South  China  subplate  is  the  Tanlu 
fault,  which  may  be  traced  from  Landsat  photographs  (Lee,  e_t  aK  1976)  northward 
across  the  border  with  the  North  China  subplate,  through  Shantung  Province,  the 
Pohai  Gulf  and  into  Liaoning  Province  (north  of  Tangshan).  The  Tangshan 
earthquake  (July  1976)  and  the  Haichang  earthquake  (February  1975)  indicate  that  the 
northern  portion  of  the  path  from  Tangshan  to  Taipei  passes  through  a  tectonically 
active  region. 

The  southern  portion  of  the  path  traverses  a  region  affected  by  the  Philippine 
Sea  plate  that  is  subducted  beneath  the  Eurasion  plate.  The  tectonics  of  this 
region  have  been  discussed  by  Sun  and  Teng  (1977).  Considerable  seismic  activity 
is  found  throughout  northern  and  eastern  Taiwan,  particularly  along  a  major  north- 
south  trending  fault  paralleling  the  east  coast  of  Taiwan.  Sun  and  Teng  (1977) 
discussed  the  importance  of  this  fault  to  a  subduction  zone  of  Mesozoic  or  Cenozoic 
age.  The  eastern  provinces  of  Fuchien,  Kiangsu,  Chiekiang,  Kwantung,  Anhwai  and 
Kiangsi  are  covered  by  large-scale  volcanic  intrusions  which  are  of  Mesozoic  age. 

These  volcanic  intrusions  are  attributed  to  a  magma  source  as  a fconsequence 
of  the  melting  from  the  subduction  of  the  Pacific  Sea  plate.  Finally,  higher  than 
normal  heat  flow  (3  to  4  heat  flow  units)  has  been  observed  north  of  Taiwan  which 
is  associated  with  areas  of  plate  subduction  (Yansuii,  et  al_.  1970).  All  these 
factors  would  contribute  to  a  shear  wave  velocity  profile  lower  than  normal  in  the 
upper  mantle. 

Conclusions 

By  using  a  multiple-filter  technique,  group  velocities  for  various  regions  of 
China  were  obtained  from  earthquakes  recorded  at  SRO  stations  in  Mashad,  Iran  and 
Taipei,  Taiwan.  For  each  sampling  path,  a  number  of  events  were  used  with  practically 
the  same  source-station  geometry  to  derive  the  dispersion  data.  An  estimate  of  un¬ 
certainty  of  the  data  was  thus  obtained.  A  nonlinear  least  square  technique  was 
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used  to  invert  the  surface  wave  group  velocity  for  shear  wave  velocity  structures. 
The  findings  concerning  the  nature  of  the  crust  and  mantle  are: 

1.  Group  velocities  obtained  for  the  pure  path  across  the  Tibetan  plateau 
are  unusually  low  when  compared  with  the  average  continental  dispersion.  Results 
of  the  inversion  indicate  that  a  4-layer  70-km  thick  crust  is  an  adequate  model 
for  the  Tibetan  plateau. 

2.  Group  velocities  for  the  mixed  path  between  Tangshan  and  Mashad,  Iran 
are  lower  than  the  average  continental  dispersion.  Results  of  the  inversion 
indicate  that  a  3-layer,  45  km  thick  crust  is  an  adequate  model. 

3.  Group  velocities  for  the  path  between  Tangshan  and  Taipei,  Taiwan  are 
closest  of  the  three  regions  to  average  continental  dispersion.  Results  of  the 
inversion  indicate  that  a  3-layer  30-km  crust  is  an  adequate  model. 

For  all  three  models,  the  structure  of  the  upper  mantle  from  information 
provided  by  the  data  is  tentative.  In  comparing  the  apriori  (SDX0)  to  the  a 
posteriori  (SOX)  standard  deviations  in  the  mantle,  it  was  observed  that  the 
final  model  was  influenced  by  the  starting  model  more  than  the  data.  Attempts 
to  increase  the  information  from  the  data  and  the  resolution  result  in  either 
unacceptably  large  standard  deviations  or  unstable  shear  velocity  models.  Only 
in  the  case  of  data  for  the  path  from  Tangshan  to  Mashad,  for  which  inversions 
were  performed  simultaneously  for  long  period  Rayleigh  and  Love  waves,  was 
adequate  information  obtained  about  the  upper  mantle. 

Previous  surface  wave  studies  of  China  have  applied  a  trial -and-error  inver¬ 
sion  approach  to  obtain  a  snear  velocity  model  from  Rayleigh  wave  group  velocity 
data.  Proper  consideration  has  not  been  given  to  either  the  information  obtained 
from  the  data,  or  the  resolution.  In  all  these  studies  the  influence  of  the 
starting  model  upon  the  final  model  was  not  discussed,  thus  the  validity  of  the 
estimated  upper  mantle  shear  velocities  presented  is  rather  questionable. 

Both  Rayleigh  and  Love  waves  should  be  inverted  simultaneously  to  obtain 
more  Information  from  the  data.  This  may  not  always  be  possible  due  to  crustal 
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anisotrophy.  Higher  mode  data  are  very  important,  but  generally  the  ex¬ 
citation  of  higher  modes  is  not  strong  for  the  shallow  earthquakes  used 
in  this  study.  In  this  regard,  future  large  earthquakes,  in  Tibet  for 
instance,  of  intermediate  focal  depths  would  be  important  sources 
for  higher  mode  data  and  they  warrant  special  attention.  Further,  the 
dispersion  curves  must  not  be  contaminated  by  surface  waves  traveling 
non-least  time  paths  or  by  complex  structure  along  the  path.  This  often 
presents  a  problem  in  tectonically  active  regions,  resulting  in  the 
inability  to  obtain  a  good  love  wave  dispersion  curve  from  the  data. 

Finally,  the  problems  due  to  epicentral  distance  must  be  considered.  If 
the  epicentral  distance  is  too  short,  proper  surface  wave  mode  separation 
will  not  be  achieved.  Only  after  those  factors  are  considered  can  one 
arrive  at  a  better  shear  wave  velocity  structure  by  means  of  surface  wave 
inversion  for  a  specific  geological  province. 
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TABLE  2 

Tangshan-Taipei  Observed  Dispersion  Data 

Rayleigh  Have  Love  Wave 

Period  (sec)  (km/sec)  (km/sec) 


64.0 

3.49 

+  .06 

3.85 

+  .15 

51.2 

3.48 

+  .09 

3.71 

+  .10 

42.7 

3.46 

+  .11 

3.69 

+  .06 

36.6 

3.48 

+  .11 

3.59 

+  .05 

32.0 

3.30 

+  .09 

3.52 

+  .08 

28.4 

3.29 

+  .05 

3.37 

+  .05 

25.6 

3.15 

+  .08 

3.31 

+  .05 

23.3 

3.04 

♦  .05 

3.19 

+  .05 

21.3 

3.00 

+  .06 

3.19 

+  .05 

19.7 

2.84 

+  .06 

3.12 

+  .09 

18.3 

2.76 

+  .09 

3.12 

+  .07 

17.1 

2.72 

+  .06 

3.04 

+  .09 

16.0 

2.67 

+  .05 

3.02 

+  .04 

15.1 

2.63 

+  .10 

2.93 

+  .05 

13.5 

2.58 

+  .13 

2.89 

+  .10 

12.2 

2.52 

+  .10 

2.85 

+  .15 

11.1 

2.59 

+  .18 

2.90 

+  .07 

10.2 

2.64 

+  .17 

2.91 

+  .15 
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TABLE  3 


Tangshan-Mashad  Observed  Dispersion  Data 


Period (sec) 


Rayleigh  Wave 
(km/sec) 


Love  Wave 
(km/sec) 


102.4 

3.73 

+  .20 

4.00 

+  .20 

85.3 

3.65 

+  .18 

3.91 

+  .13 

73.1 

3.69 

+  .10 

3.91 

+  .10 

64.0 

3.57 

+  .10 

3.85 

+  .05 

56.9 

3.65 

+  .18 

3.74 

+  .06 

51.2 

3.54 

+  .05 

3.70 

+  .11 

46.5 

3.47 

T  .05 

3.66 

+  .10 

42.6 

3.43 

+  .05 

3.56 

+  .06 

39.4 

3.35 

+  .05 

3.46 

+  .05 

36.6 

3.22 

+  .05 

3.38 

+  .05 

34.1 

3.14 

+  .04 

3.36 

+  .05 

32.0 

3.10 

+  .05 

3.31 

+  .05 

30.1 

3.01 

+  .06 

3.29 

+  .06 

26.9 

2.93 

+  .06 

3.21 

+  .05 

24.4 

2.91 

+  .08 

3.18 

+  .05 

22.3 

2.85 

+  .10 

3.16 

+  .07 

20.5 

2.74 

+  .06 

3.17 

+  .05 

18.9 

2.78 

+  .10 

3.14 

+  .05 

16.0 

2.82 

+  .18 

3.06 

+  .06 

13.8 

2.75 

+  .07 

3  -  0 1  ?+  .06 

12.2 

2.79 

+  .08 

2.95 

+  .10 

10.9 

2.69 

+  .08 

2.89 

±  *15 
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TABLE  4 


Yunnan-Mashad  Observed  Dispersion  Data 


Rayleigh  Wave (km/sec) 

Period (sec)  _ Fundamental  First  Higher  Mode 


102.4 

3.52  +  .15 

85.3 

3.49  +  .15 

73.1 

3.41  +  .10 

64.0 

3.25  +  .15 

56.9 

3.12  +  .10 

51.2 

3.05  +  .10 

46.5 

3.00  +  .10 

42.7 

2.95  +  .10 

39.4 

2.87  +  .07 

36.6 

2.79  +  .07 

34.1 

2.73  +  .07 

32.0 

2.71  +  .07 

30.1 

2.73  +  .15 

27.0 

2.67  +  .15 

24.4 

2.85  +  .15 

22.3 

2.79  +  .05 

20.5 

2.79  +  .05 

4.09  +  .15 

18.9 

2.77  +  .05 

3.80  +  .17 

16.0 

2.77fc+  .06 

3,-48f+  .10 

13.8 

2.77  +  .06 

12.2 

2.64  +  .15 

3.17  +  .08 

10.9 

2.71  +  .10 

3.15  +  .07 
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TABLE  S 


Szechwan-Mashad  Observed  Dispersion  Data 
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Rayleigh  Wave  (Jcm/sec) 

_ Fundamental  First  Higher  Mode 
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TABLE  8  :  Model  Fit  for  Path  Yunnan-Mashad 


TABLE  10:  Model  Fit  for  Path  Tangshan-Mashad 


U  i  Resolution  Matrix  for  Path  Tan 


TABLE  12 :  Model  Fit  for  Path  Tangshan-Taipei 


TABLE  13:  Resolution  Matrix  for  Path  Tangshan-Taipei 


Figure  1.  Regional  tectonic  subdivision  of  China  with  inferred  crustal  thickness  shown  in  5  km 
contours  (after  Chinese  Academy  of  Sciences,  1974). 
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Figure  2.  Surface-wave  paths  (triangles  are  earthquakes,  dots  are  SRO  stations) 


Figure  3.  Unrotated  seismograms  for  the  path  Tangshan  -  Mashad, 
August  8,  1976  (from  top  to  bottom:  vertical,  NS,  and 
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Figure  5.  Results  of  multiple  filtering  of  Rayleigh  wave  vertical  component,  Tangshan  -  Mashad, 
August  8,  1976. 
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Figure  7.  Results  of  multiple  filtering  of  Love  wave,  Tangshan  -  Mashad,  August  8,  1976. 
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Figure  10.  Rayleigh  wave  group  velocity  for  three  events  of  the 
path,  Yunan  -  Mashad. 
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Figure  11.  Group  velocity  of  Rayleigh  waves  for  five  events  of  the 
path,  Tangshan  -  Mashad. 
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Figure  13.  Vertical  component  Rayleigh  wave  group  velocity  for  four 
events  of  the  path,  Tangshan  -  Taipei. 


Figure  14.  Radial  component  Rayleigh  wave  group  velocity  for  three 
events  of  the  path,  Tangshan  -  Taipei. 
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Comparison  of  Rayleigh  wave  group  velocities  for  several  paths  across  China  with 
average  continental  dispersion  (Ewing  et  aL »  1957). 
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Figure  19.  Resulting  shear  wave  velocity  model  for  the  parth  Szechwan-Mashad 
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Figure  21.  Comparison  of  the  Canadian  Shield  model  (Brune  and  Dorman,  1963)  and  the 
Gutenberg  model  (Takeuchi  et  aK »  1964)  with  model  SMI  ana  YMl . 
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Figure  22.  Theoretical  and  observed  group  velocities  for  the  Tangshan-Mashad  path. 
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Figure  23.  Resulting  shear  wave  velocity  model  for  the  Tangshan-Mashad  path. 
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Figure  26.  Resulting  shear  wave  velocity  for  the  path  Tangshan-Taipei . 


Love  Wave 


Observed  and  theoretical  group  velocities  for  the  path  Tangshan-Taipei 
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Comparison  of  Canadian  Shield  Model  (Brune  and  Dorman,  1963)  and  the 
Gutenberg  Model  (Takeuchi  et  al.,  1964)  with  model  TT1. 


